Effective mass enhancement in two-dimensional electron systems: the role of 

interaction and disorder effects 
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Recent experiments on two-dimensional (2D) electron systems have found a sharp increase in the 
effective mass of electrons with decreasing electron density. In an effort to understand this behavior 
we employ the many-body theory to calculate the quasiparticle effective mass in 2D electron systems. 
Because the low density regime is explored in the experiments we use the GWT approximation where 
the vertex correction F describes the correlation effects to calculate the self-energy from which the 
effective mass is obtained. We find that the quasiparticle effective mass shows a sharp increase with 
decreasing electron density. Disorder effects due to charged impurity scattering plays a crucial role 
in density dependence of effective mass. 

PACS numbers: 71.10.-w, 71.10. Ca, 71.30.+h 
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There has been a large amount of experimental and 
theoretical activity in recent years to understand the 
ground state properties of homogeneous two-dimensional 
(2D) electron systems. Advances in fabrication tech- 
niques have made it possible to probe various quanti- 
ties of interest in high quality and very low density sam- 
ples. Most notably, the observation of a metal-insulator 
transition[l| in these systems provides a major motiva- 
tion to study the various physical properties. In re- 
cent experiments the spin susceptibility, Lande g-factor, 
and effective mass are measured for 2D electron sys- 
tems made of Si-MOSFETS and GaAs quantum-well 
structures. 

H H E3 IE II In particular, Shashkin et 
al. 0, 3 reported a sharp increase of effective mass near 
the critical density at which the system starts to show de- 
viations from the metallic behavior. On the other hand, 
Pudalov et al. 0] have found only moderate enhancement 
of the spin susceptibility and effective mass in their sam- 
ples. 

There has been a number of calculations of the quasi- 
particle properties including the effective mass of 2D elec- 
tron gas employing a variety of approximations. 8] More 
recent theoretical calculations of the effective mass of 2D 
electrons concentrated on the density, spin polarization 
and temperature dependence. 0, [Hi HU, H3, UM 

In view of the different experimental results and their 
controversial interpretation we have addressed in this 
work the density dependence of the effective mass in 
an interacting electron system at T = in the pres- 
ence of charged impurities. We employ the GWT 
approximation 14) which includes the vertex corrections 
in an approximate way to calculate the quasiparticle ef- 
fective mass. The local-field factor describing the correla- 



tion effects which enters the vertex function T is obtained 
within the memory function formalism and the self- 
consistent field method. We had previously shown [Fo| 
that such an approach correctly describes the anoma- 
lous behavior of the thermodynamic compressibility of 
2D electron systems as reported experimentally, fig In 
the present work, we extend our earlier considerations to 
calculate the effective mass. We find that the quasipar- 
ticle effective mass is greatly enhanced at low density in 
the same region when the compressibility diverges. 

In the following we first outline the theoretical frame- 
work with which we calculate the quasiparticle effective 
mass of 2D electron system. We next present our results 
within various levels of approximations. We discuss the 
results of our calculations in view of other theoretical ap- 
proaches and experimental findings. We conclude with a 
brief summary. 

We consider a 2D electron system interacting via the 
long range Coulomb interaction V q — 2ire 2 /(eoq) where 
e is the background dielectric constant. The system is 
characterized by the dimensionless interaction strength 
r s = 1/ 1 (■Kna* B 2 ) 1 / 2 , where n is the 2D electron density 
and a* B = H 2 eo/(me 2 ) is the effective Bohr radius defined 
in terms of the band mass m of electrons in the semicon- 
ductor structure. 

We use the theoretical framework developed by Thakur 
et a/.|l7| employing the memory- function formalism and 
the self-consistent field method to calculate the density- 
density response function of a disordered electron sys- 
tem. The effect of disorder is to dampen the charge- 
density fluctuations and modify the response function. In 
a number-conserving approximation the density-density 
response function for noninteracting electrons is given by 
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where 7 is the scattering rate. The correlation effects 
are described by the generalized random-phase approxi- 
mation (RPA) for the interacting system density-density 
correlation function 
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in which the static local-field factor G(q) embodies the 
correlation effects. In this work we use the self-consistent 
field method of Singwi et to calculate G{q). As a 

simplified model we also con sider the Hubbard local-field 
factor given as G(q) = q/2y / q 2 + k 2 F . 

Within the memory function formalism the scattering 
rate is expressed in terms of the screened disorder poten- 
tial and the relaxation function asj^J IHJ 
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where Ui mp (q) — V q e~ qd is the impurity potential for 
charged impurities located at a distance d away from the 
2D electron layer, and rij is the impurity concentration. 
The relaxation function is given bv[rH IT^| </>o(<M7) = 
[X(ff,*7) - Xo(ff)]/(»7) and e(q) = 1 - V q [l - G(q)} Xo (q) 
is the static screening function. Because the scattering 
rate 7 depends on the screening function e(q) which itself 
is determined by the disorder included response function 
the above set of equations are solved self-consistently. 

The quasiparticle properties of the 2D electron system 
are obtained from the self-energy function[l!] E(fe,w) 
which we calculate at zero temperature. Since we are 
interested in exploring the interaction effects we include 
the vertex corrections to the self-energy and employ the 
GWT approximation.^) The self-energy in the GWT 
approximation is written as a sum of two terms £(fc, u>) — 
£ line 0, uj) + £ polc (fc, uj) where 
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in which = k 2 /2m — Ef is the single-particle energy 
measured relative to the Fermi energy. The vertex func- 
tion in the local-approximation is given as|l4| 
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in terms of the local-field factor G(q) describing cor- 
relation effects beyond the RPA. The dielectric func- 
tion appearing in Eqs. (4) and (5) is given by e(q,ui) = 
1 — V q x (q, uj; 7)r(g, lu). The above expressions for the 
self-energy reduce to the GM^-RPA results when we set 
G(q) = 0. Furthermore, taking 7 = 0, we recover the 
results for a clean system. 

We use the on-shell approximation to the self-energy in 
the single particle spectrum = + S(fc,£fc) to obtain 
the effective mass perturbativelv|8l. Ull H4| 
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where the frequency and momentum derivatives of 
E(fc,w) are evaluated at the Fermi surface. Such an ap- 
proach is argued to be more appropriate over solving the 
full Dyson's equation Eh = £fc + S(fc, Ef.), since the self- 
energy is calculated using the noninteracting Green func- 
tion. The resulting scheme incorporates the higher order 
diagram contributions better. |8|, LUl LL|| 

In the numerical calculations we specialize to GaAs 
systems for which some measurements of the effective 
mass are undertaken. Since the dominant scattering 
mechanism is known to be that due to the charged im- 
purities we take d = 250 A for the setback distance 
in Ui mp (q), and consider rii to be of the order of ~ 
10 10 cm -2 . We first solve the self-consistent equations for 
the scattering rate 7 and the local-field factor G(q) as a 
function of r s and impurity density m. These quantities 
determine the dynamic screening function e(q, w) and the 
vertex function T(q, ui). Afterward the self-energy and its 
derivatives are evaluated to find the effective mass m* jra. 

In Fig. 1 we show the effective mass m*/m as a func- 
tion of r s calculated in various theoretical approaches. 
We observe that even at these relatively higher densities 
there are notable differences between the GW and GWT 
approximations indicating the importance of vertex cor- 
rections. It appears that the correlation effects suppress 
the effective mass renormalization. The disorder effects 
due to charged impurity scattering tend to increase the 
effective mass with respect to the clean system. There is 
still a significant reduction in m*/m, however, compared 
to the GI^-RPA result. 

We next display the quasiparticle effective mass at 
much lower densities. Strictly speaking, GW approxima- 
tion is valid only in the high density limit, inclusion of 
the vertex corrections in GWT approximation improves 
the regime of validity. In any case, we wish to explore the 
effective mass trends in the low density, strongly interact- 
ing region. Figure 2 shows m*/m calculated in various 
theoretical approaches. GW-KPA yields a modest en- 
hancement for the whole density range. The suppression 
found in Fig. 1 for the GWT approximation reverses its 
behavior around r s rj 4 and shows an enhancement rel- 
ative to the GW^-RPA results at lower densities. Qual- 
itatively similar behavior is obtained when we use the 
simple Hubbard local-field factor within the GHT ap- 
proximation. A notable feature of the GWT approxi- 
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FIG. 1: The quasiparticle effective mass as a function of r s in 
the range < r s < 2. The solid and dot-dashed lines indicate 
the GW-KPA and GWT approximations, respectively. GWT 
approximation which uses the Hubbard local-field factor is 
indicated by the dotted line. GWT approximation including 
the charged impurity scattering (with impurity concentration 
rii = 0.5 x 10 10 cm -2 ) is shown by the dashed line. 
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FIG. 2: The quasiparticle effective mass as a function of r s in 
the range < r s < 10. The solid and dot-dashed lines indicate 
the GW-KPA and GWT approximations, respectively. GWT 
approximation which uses the Hubbard local-field factor is 
indicated by the dotted line. GWT approximation including 
the charged impurity scattering (with impurity concentration 
rii = 0.5 x 10 10 cm -2 ) is shown by the dashed line. 



mation results is that effective mass exhibits a sharp in- 
crease around r s ~ 8. That the strong interaction effects 
would lead to a large enhancement in to* /to is also evi- 
dent when the Hubbard local-field factor is used within 
the GWT approximation. Finally, when charged impu- 
rity scattering effects are included in the calculation we 
find that a similar sharp increase in to* /to occurs at a 
smaller r s value. We have also calculated the effect of im- 
purity scattering for different parameter values of d and 
rii and found qualitatively similar results. 

Although the results of GW and GWT approximations 
at large r s are should be taken in with caution, the low 
density trends of to*/to should be indicative. In this per- 
spective our calculations indicate that the effective mass 
enhancement in 2D electron systems can be accommo- 
dated within the Fermi liquid theory when the vertex cor- 
rections describing the strong correlation effects are taken 
into account. In particular, a sharp increase in to*/to as 
shown in Fig. 2 is quite suggestive in view of the recent 
experimental findings. U Effective mass enhancement 
is also observed in 2D neutral Fermi systems pfl| and rep- 
resented by a GW -type calculation. [2l| Our calculations 
also show that the r s value at which to* /to exhibits a 
sharp increase can be controlled by disorder effects. In a 
self-consistent scheme where remote charged impurities 
are taken into account we find that the large enhance- 
ment in to* /to occurs at a higher density compared to 
the strongly interacting clean system. 

Recent theoretical approaches [ij 0] have modeled the 
low density electron liquid as close to the Wigner crys- 
tallization to obtain a strong increase in the effective 
mass. On the other hand, Morawetz found a di- 
vergent behavior in to* /to at the metal-insulator tran- 
sition by considering the scattering from heavy impurity 
ions, and Galitski and Khodel^| attribute divergence of 
the effective mass to the density wave instability. Our 
theoretical scheme considers the metallic regime, there- 
fore we cannot distinguish the nature of the possible new 
state beyond the critical r s value. However, our earlier 
calculations^^ of the anomalous behavior of the com- 
pressibility at around the same range of r s values points 
to a possible metal-insulator transition approaching from 
the metallic side. A related quantity of interest would be 
the spin susceptibility or the g factor for which exper- 
imental results are available. These quantities require 
the calculation of density-density response function and 
the local-field factor as functions of the spin polarization 
which were not undertaken in this work. 

Our calculations were performed at T = and for zero 
thickness 2D electron layers. It would be interesting to 
extend our work to finite temperatures and to finite width 
quantum wells to make better contact with experiments. 
As the Coulomb interaction effects will be less strong in 
quantum wells, it is expected that the enhancement of 
to* /to will be less marked. 

In summary, within a many-body approach which 
takes the electron-electron and electron-impurity in- 
teraction effects into account we have calculated the 
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effective mass of a 2D electron system at zero tempera- 
ture. We have found within the commonly used on-shell 
approximation that m* /m is highly enhanced at larger 
values of the density parameter r s as a result mainly 
of the correlation effects. The interplay between the 
correlation effects and the impurity scattering influences 
qualitative changes in this behavior. Our comparative 



results should provide insight into the workings of 
many-body methods in the strong interaction regime. 
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